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We present ac susceptibility and specific heat measurements taken on samples of LiHo I Yi_ a: F4 in 
the dilute limit: x — 0.018, 0.045, 0.080 and 0.12. Susceptibility measurements show glassy behavior 
including wide absorption spectra that continually broaden with decreasing temperature. Dynamical 
scaling analyses show evidence of finite-temperature spin-glass transitions, the temperatures of which 
match those of recent theoretical work. A surprisingly long intrinsic time constant is observed in 
these samples and is found to be inversely correlated with the concentration of magnetic moments, x. 
Our results support the picture that this behavior is largely a single-ion effect, related to the random 
transverse fields generated by the off-diagonal component of the dipolar interaction and significantly 
slowed by the important nuclear hyperfine interaction. Specific heat measurements show broad 
features due to the electronic spins on top of a large Schottky-like nuclear contribution. Unusually, 
the peak position of the electronic component is found to be largely concentration independent, 
unlike the glass transition temperature. 



I. INTRODUCTION 

The compound LiHoF4 and derivatives obtained 
through dilution of the magnetic Ho 3+ moments with 
non-magnetic Y 3+ , have long been considered to be excel- 
lent representations of the dipolar Ising mo del and have 
been studied extensively as model magnetsP^ In recent 
years, however, this series of materials has become the 
subject of considerable debate, both experimental^ and 
theoreticalpESI particularly regarding the dilute (small 
x) limit and the question of whether the introduction of 
randomness and frustration originating from the dipolar 
interaction leads to a spin-glass transition or something 
more exotic EEH 

There are two main components to this debate. First 
is the largely theoretical question of whether the sim- 
ple toy model that is thought to describe LiHo x Yi_ a; F4, 
the dilute, dipolar-coupled Ising model, ought to ex- 
hibit a finite-temperature spin-glass transition or not. 
While mean field theory suggests that there is a spin- 
glass transition, 13 some Monte Carl o sim ulations have 
called that conclusion into question ! 8 * 9 * 14 * However, re- 
cent Monte Carlo results, implementing parallel temper- 
ing have shown strong evidence of a finite-temperature 
spin-glass transition.^ 

The other major part of the debate deals with the real 
material and the possibility that the toy model is in fact 
not an adequate description of the system's underlying 
model. Several experimental groups have come to dif- 
ferent conclusions regarding the dilute regime. One re- 
search group has come to the conclusion that at x = 0.167 
and x — 0.20, the materials undergo a spin-glass transi- 
tion at 130 mK and 150 mK, respectively, as determined 
from the temperature dependence of the ac su sceptibility 
x(w) and nonlinear susceptibility ygl 1 * 6 * 15 !^ whereas at 
x = 0.045 an unusual spin liquid state appears? 1 * 11 * 12 * 17 *^ 
This spin liquid or "antiglass" state is primarily charac- 
terized by a narrowing of the absorpt ion spectrum x"^) 
as the temperature is lowered!- 1 * 11 * 17 * Several other inter- 



esting effects were observed in the 4.5% material includ- 
ing shar p features in the specific heat at 120 and 300 
m |£jUii| contrasting with a smooth T -0 - 75 power law de- 
pendence of the dc susceptibility,^ hole burning in the 
absorption spectrum and possible coherent oscillations.^ 
A second experimental group, in contrast, has found 
qualitat ivel y similar beh avio r between 4.5% and 16.7% 
samples.^! Jonsson et alW^ have concluded that neither 
sample shows a sufficiently diverging \3 to be considered 
a spin glass with a finite-temperature glass transition, 
T g , instead suggesting that the materials are a kind of 
superparamagnet with thermally activated dynamics. 

We present here, and in other short articles'^ a third 
point of view: that there is a finite-temperature spin- 
glass transition in samples of LiHo x Yi_ :E F4 for x < 0.2, 
and that there is no exotic antiglass state at x = 0.045. 
We have performed ac susceptibility measurements us- 
ing a specially designed SQUID magnetometer on sam- 
ples of x — 0.018, 0.045 and 0.080 and have performed 
specific heat measurements on those same samples plus 
an additional x = 0.12 sample. The data present sev- 
eral puzzles and conclusions are hampered in this system 
by extremely long time constants, even well above T g , 
which can be seen to be a result of single ion physic s and 
the importance of the nuclear hyperfine interaction! 19 * 20 * 
Nonetheless, the evidence largely supports the picture of 
spin glass physics in these materials. 

This article begins with a discussion of the micro- 
scopic Hamiltonian that is thought to describe the 
LiHo 2; Yi_ 2 ;F4 system and a review of past work on the 
phase diagram of this series of materials. We concen- 
trate on the low-x part of the phase diagram, but refer 
the reader to the recent and more inclusive review article 
Ref.[2T] This is followed by a description of the apparatus 
and methods used in this work, specifically, our specific 
heat and ac susceptibility experiments. In this work, we 
expand on previous measurements of the specific healP 
and ac susceptibility^ of LiHo 2: Yi_ 2 ;F4, covering a larger 
range of x. Specifically, the new work contained within 
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operators,^ as 
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FIG. 1: (color online) The crystal structure of L1H0F4. 
The Ho 3+ moments may also be randomly replaced by non- 
magnetic Y 3+ ions during crystal growth. Only the F~ ions 
immediately surrounding one of Ho 3+ /Y 3+ sites have been 
drawn in order to simplify the image. The crystal structure 
provides a crystal field with S4 symmetry around the Ho 3+ 
ions, leading to an Ising ground-state doublet. On dilution, 
random frustration is provided by the anisotropy of the dipo- 
lar interaction. The dipolar interaction is ferromagnetic be- 
tween nearest neighbors (NN) but antiferromagnetic between 
next nearest neighbors (NNN). 



this article consists of ac susceptibility measurements on 
x = 0.018 and x — 0.080 samples as well as a specific 
heat measurement at x — 0.12. In Section IV, we present 
our results along with scaling analyses that are stron gly 
suggestive of spin glass physics. An explanatio n 1 19 ! 20 ! for 
some of the dynamical behavior that is observed, based 
on the single-ion Hamiltonian of the system, is presented 
in Section V. Our results are then compared with those 
of other groups and disparities in data and interpretation 
are discussed. 



II. BACKGROUND 

The magnetic Ho 3+ ions in LiHo x Yi_ l: F4 possess spin 
5 = 2, orbital angular momentum L — 6 and total an- 
gular momentum J — 8 moments as a result of Hund's 
rules and strong spin-orbit coupling, typical of rare earth 
ions. The J = 8 levels are further split in energy by 
the crystalline electric field, of S4 symmetry. The result- 
ing Hamiltonian can be expressed in terms of Steven's 
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Ho 3+ is a non-Kramers ion, so it is only by virtue of 
the symmetry of the crystal field, that an Ising doublet 
ground state results. This ground-state doublet, consist- 
ing of states I f) and | i) , has no matrix elements of J x 
or J y . It does, however, exhibit large moments along the 
z-direction (c-axis). 

(tlJ.lt) 
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thus the effective g-factor is 
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where g,j = 5/4 is the Ho 3+ Lande g-factor. The next 
excited state of the Ho 3+ ions, I7), has been measured 
to be roughly 11 K in en ergy a bove the Ising doublet by 
several groups in the past 122121] Thus for the temperatures 
of interest, say near and below T c (x = 1) = 1.53 KpmUD 
this system is very well described as an Ising magnet. 

Because the 4/-electrons of the rare-earth ions are 
tightly bound, they tend to possess relatively small 
exchange interactions. When this is combined with 
the large magnetic moments, the dipole-dipole inter- 
action can become energetically quite important. In 
LiHOa-Y^^F^ the dipolar interaction is in fact domi- 
nant; the nearest neighbor exchange interaction has been 
estimated experimentally to be roughly half the nearest 
neighbor dipolar interaction^ and in a recent theoretical 
work was determined to be a factor of 10 smallerP The 
dipole-dipole Hamiltonian is given as 
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Within the Ising doublet manifold, however, it can be 
written instead as 



nj \ ~> Mo 2 2 
Ud = 1^ 4^effMB 
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where the Sf are spin- 1/2 operators and 
g cS = 2gj(t\J z \t)- 



(5) 
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This interaction has two important distinctions. First, 
it is long range in nature, falling off as 1/r 3 . This means 
that every spin is coupled to every other spin, at least 
to some degree, making a percolation threshold impos- 
sible, for example. Second, the sign of the interaction 
is dependent on the angle of the vector connecting the 
spins and this always leads to a degree of frustration 
in three dimensions. If the vector connecting spins is 
aligned with the Ising c-axis (so that they are stacked on 
top of each other in Fig. [T]), they are ferromagnetically 
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coupled; if they are positioned within the same a6-plane 
(perpendicular to the Ising axis), they are antiferromag- 
netically coupled. In LiHo K Y 1 _ x F4, nearest neighbors 
(NN) exhibit ferromagnetic coupling where next near- 
est neighbors (NNN) exhibit antiferromagnetic coupling, 
as shown in Fig. [IJ causing frustration. The ferromag- 
netic couplings dominate in the pure material, which has 
a ferromagnetic transition temperature of 1.53 K!212itf!H 
Upon randomly diluting the system with non-magnetic 
Y 3+ , however, this frustration is exposed and is one of 
the main ingredients that gives rise to the glassy physics 
and possible spin glass state seen at x < 0.25. 1 

Introduction of a magnetic field results in the Zee- 
man energy Hz — ffJMfiH • J. Much of the past re- 
search performed on this system has aimed at under- 
standing the effect of quantum fluctuations, introduced 
with a magnetic field transverse to the c-axis, or H±. 
At the single-ion level, such a transverse field splits the 
Ising doublet by an energy of A(H±) — 2T and mixes it 
with the next excite d state I7). When projected onto a 
5=1/2 model j22HHD one has an effective transverse mag- 
netic field of r so the contribution TS X . The result is a 
manifestation of the famous transverse field Ising model 
(TFIM)Pl which is one of the simplest models known 
to exhibit a zero-temperature, quantum phase transition 
(QPT). This paradigm was applied by Bitko et al. to 
the ferromagnetic parent compound {x = 1)P As H± is 
increased, quantum fluctuations become more and more 
powerful and eventually melt the ferromagnetic order, 
leading to a quantum paramagnetic state. In LiHoF 4 , 
this is seen to occur at a critical field Hi = 4.9 TP A 
number of theoretical works have attempted to fit the 
(H± , T) phase boundary, but a p ersistent quantitative 
disagreement has been foundPSESI The previous experi- 
mental phase boundary has recently been confirmed with 
dilatometry measurements,^ and the mismatch between 
theory and experiment remains an outstanding problem. 

The last component to the system's Hamiltonian to 
mention is the hyperfine coupling to the / = 7/2 nuclear 
moment. Holmium, in fact, has an extraordinarily large 
hyperfine interaction described by 

•Hhf = AL J (7) 

where A = 40.21 mK! 27 | 35 | 3 ^l In the perfect Ising model, 
where no excited crystal field states are considered, this 
leads to 8 distinct electronuclear energy levels, separated 
by 207 mK. While smaller than 1.53 K, thus not terri- 
bly influential in zero field for the parent compound, it is 
easy to see that the hyperfine coupling can become ener- 
getically important as one moves to low concentrations of 
x < 0.25, where glassy physics occurs. The hyperfine in- 
teraction also has important consequences near the QPT 
of higher concentration samples, creating a non-trivial 
, T) phas e diag rarrP* 3 ^ and influencing the spectrum 
of excitations) 23 1 37 ' 

While a number of very interesting experiments have 
been performed on the ferro magn etic stoichiometries x = 
0.46, x = 0.67 and x = lJUSHU] we concentrate here 



on the dilute glassy portion of the phase diagram at 
x = 0.2 and below. The first experiments performed 
on an x — 0.167 sample 1 showed largely conventional 
spin glass behavior. Glassy relaxation was observed in 
the ac susceptibility, x(/), and a broad feature was ob- 
served in the specific heat C . Most evidence seemed to 
point to the existence of a dipolar spin-glass transition as 
was anticipated theoretically by Stephen and AharonyP^ 
LiHo a; Yi_ 2 ;F4 is not the first system to be investigated 
as a possible dipolar spin glass. Dilute EuxSri-^S was 
also considered as a candidate for such behavior, but is 
complicated by strong exchange interactions.^ The low x 
range of the LiHo ;r Yi_. E F4 series represents an important 
test of theories that either propose or deny the existence 
of a finite glass transition in a purely dipolar system. 

The effects of an applied transverse magnetic field were 
also tested at the x = 0.167 stoichiometry with the goal 
of studying a quantum phase tra nsitio n from spin glass 
order to quantum paramagnetismP^^ This was achieved 
using linear ac susceptibility^ and also nonlinear suscep- 
tibility, %3, measurements™ The results were, seemingly, 
an (H± , T) phase boundary with a classical phase tran- 
sition at T g ~ 130 mK and a quantum phase transition 
at He — 1.2 T. Although at higher temperatures (lower 
fields), X3 is found to exhibit a very sharp peak at the 
phase transition, at lower temperatures (higher fields), 
the X3 feature become s extremely broadened. Recent 
theoretical worlPSMElHlll nas come to explain this effect 
through static random fields that are induced by the ap- 
plied transverse field. In fact, for any finite H±, there 
is no longer a true spin glass state, but only domains 
of spin glass order that decrease in size with transverse 
field. The phase boundary is, thus, a crossover rather 
than a true line of phase transitions. However, for low 
fields where it is quite sharp, this distinction is not easily 
made experimentally. 

At lower concentration, x = 0.045, Reich et a/pmi 
observed an unusual narrowing of the absorption spec- 
tra x"(w) with decreasing temperat ure. This is at odds 
with typical behavior of spin glasses^EH a nd, in partic- 
ular, with the clear broadening of x"^) that is seen at 
x = 0.167. 1 -' This stoichiometry has thus been referred 
to as an "antiglass" . In later work by Ghosh et al. , the 
same research grou p, significantly more unusual results 
were observed! 11 * 12 ' In those results, the absorption spec- 
tra narrow appreciably at low temperatures and develop 
a strong asymmetry 11 Hole-burning experiments were 
also performed where the system was saturated in the 
vicinity of a pump frequency, resulting in a hole in the 
spectrum as obtained with a variable probe frequency.^ 
Most surprisingly, after cutting off the oscillating applied 
field, a ringing of the magnetization was observed, decay- 
ing over several seconds, reminiscent of coherent oscilla- 
tions or a free induction decayP^ Additionally, instead of 
a single broad feature in the specific heat, two sharp fea- 
tures were observed at around 100 mK and 300 mK and 
a very small percentage of the expected R In 2 entropy is 
accounted for with a numerical integral of C/TP21 This 



4 



0.018 



0.045 



0.080 



■o 0.08 
O 




10* 10 

/(Hz) 



10"' 10" 10' 

/(Hz) 



10 10 

/(Hz) 



FIG. 2: (color online) Frequency scans of the ac susceptibility of LiHo I Yi_ I F4 for x — 0.018, x — 0.045 and x — 0.08. The 
top panels show the in-phase susceptibility or x'if) while the bottom panels show the absorption spectra x"(f)- Temperature 
increases from left (blue) to right (red). These data have been corrected for the demagnetization effect. 



unusual C(T) and a T -0 75 power law seen in the dc limit 
of x were proposed to result from quantum entanglement 
of pairs of Ho 3+ moments.^ 

In more recent years, the glassy portion of the phase 
diagram has sparked a great deal of controversy. Re- 
cent specific heat measurements^ and ac susceptibility 
results^, by the authors of this work, at x = 0.045, have 
not found any of the exotic antiglass features that were 
previously reported. Instead, results consistent with spin 
glass physics were observed^. These results and others 
will be discussed in Section [TV] 

Jonsson et alW measured both an x = 0.045 sample 
and an x = 0.167 sample, performing linear and non- 
linear susceptibility measurements. While not finding 
any major qualitative differences between those two sam- 
ples, thus not observing a dramatic change to an antiglass 
state, they also argue that there is not sufficient evidence 
for a spin glass state, in either of the samples. This find- 
ing opened up a debate regarding the size of magnetic 
field and sweep rate with which X3 should be measured 
in order to make accurate concl usio ns regarding the ex- 
istence of a spin-glass transition!^ 

From a theoretical point of view, it seems that the de- 
bate is largely settledPil2!J The Monte Carlo simulations 
that did not reveal a spin-glass transition 8 9 focused on 
the Binder ratio g as the mark of spin glass freezing. 
However, Tarn and Gingraa^ found that the correlation 
length £sq was a much better probe of a freezing transi- 
tion as in conventional Edwards- Anderson spin glasses.^ 
Through a finite size scaling analysis, they showed a 
probable existence of a finite T g in a model largely rep- 
resentative of LiHo a; Y 1 „ ;I .F4. Schechter and StamppSl 
meanwhile, have discussed the effects of off-diagonal com- 



ponents of the dipolar interaction and the nuclear hyper- 
fine coupling and their introduction of quantum fluctua- 
tions, and have concluded that they are not sufficient to 
suppress spin glass ordering. 



III. EXPERIMENT 

High quality single crystals of LiHo I Yi_ a; F4 were ob- 
tained commercially^ and their characterization has 
been discussed previously.^ Different parts of the same 
single crystals were used in Refs. 13151511521 Here, we have 
performed ac susceptibility measurements oni = 0.018, 
x = 0.045 and x — 0.080 samples, whereas we present 
specific heat measurements on four different samples: 
x = 0.018, 0.045, 0.080 and 0.12, the first three of which 
were measured previously in Ref . |3] 

The ac susceptibility measurements described here 
were performed on a SQUID magnetometer, designed 
and implemented with glassy magnetic materials in mind, 
chosen for its flat response over many decades of fre- 
quency. Where conventional susceptometers depend on 
inductive pick-up and the signal is proportional to fre- 
quency, a SQUID magnetometer, with a superconduct- 
ing flux transformer, uses the Meissner effect which is 
frequency independent. Additionally, the extremely high 
magnetic field sensitivity of a SQUID permits the use of 
excitation and pickup coils with very few turns, thereby 
avoiding resonances until frequencies well into the MHz 
range. Since only small magnitude excitations need to be 
employed, significant heating within the sample can be 
fairly easily avoided. A less than 20 mOe peak-to-peak 
oscillating magnetic field was employed in these measure- 
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ments. 

Our particular magnetometer design consists of a 375- 
turn NbTi primary coil and a second-order (5-7-5 turn) 
gradiometer as a pickup coil. The gradiometer is made 
from Nb wire and is connected to a hand- wound Nb input 
coil to the SQUID via superconducting contacts, forming 
a flux transformer. A NbTi trim coil (in parallel with the 
excitation coil) is included at one end of the gradiometer 
to adjust the balance of the circuit. Coils are wrapped 
on phenolic coil forms which are press fit together con- 
centrically. 

Several layers of noise shielding are employed. Pb su- 
perconducting shields surround the SQUID and the en- 
tire magnetometer, respectively. The 1 K radiation shield 
of the cryostat is plated with Pb, and two /i-metal shields, 
one in the He bath and one outside the dewar, are em- 
ployed. The /i-metal shields also ensure that the sample 
and SQUID are cooled in extremely low fields (expected 
to be less than 5 mOe). The SQUID is run in a flux 
locked loop with a 100 kHz modulation signal.^ 

Samples are heat sunk to a sapphire rod which is in 
turn heat sunk to the mixing chamber of a 3 He/ 4 He di- 
lution refrigerator. Most of the susceptibility data pre- 
sented here were taken on samples cut to be needle- 
shaped, with the Ising or c-axis along the length of the 
crystal. Specifically, the needle-like sample geometries for 
the 1.8%, 4.5% and 8.0% samples were 0.61 x 0.89 x 5.8 
mm 3 , 0.57 x 0.77 x 7.7 mm 3 and 0.66 x 0.94 x 7.44 mm 3 , 
respectively. Because of the large magnetic moments on 
the Ho 3+ sites, the demagnetization effect is very impor- 
tant to consider in these materials, particularly as x gets 
larger. Specifically, the demagnetization effect mixes the 
the in-phase and out-of-phase susceptibilities, resulting 
in an appreciable phase shift at lower temperatures and 
distorting the measured temperature dependence of the 
time constants of relaxation in the system.^ This makes 
it crucial to obtain a good calibration of the susceptibility 
and then correct for demagnetization. 

For the 4.5% and 8.0% samples, calibration was ob- 
tained by comparing with different sample geometries 
(0.57x0.77x3.3 and 0.66x0.94x2.6 mm 3 respectively). 
An overall amplitude calibration pre-factor for each ge- 
ometry was adjusted so that the results overlapped well 
after the demagnetization correction was applied. For a 
more in-depth description of this method of demagneti- 
zation correction, see the appendix of Ref. 1551 In the 
case of the 1.8% sample, where the moments are suffi- 
ciently dilute that the demagnetization correction is not 
significant, calibration was performed by comparing to 
a superconducting Pb calibration standard. A third ge- 
ometry of the 4.5% sample (0.57 x 0.77 x 1.2 mm 3 ) was 
also measured in order to thoroughly test for possible ge- 
ometry effects, beyond those expected from the standard 
demagnetization effect. 

Specific heat measurements, described in some detail 
previously^, were performed using the quasi-adiabatic 
heat pulse technique with a long time constant of equi- 
libration. A Ru02 thermometer and metal-film heater 



are fixed directly to the samples, avoiding the use of a 
substrate and possible heat leaks. Leads to the heater 
and thermometer are made from 6 fim diameter, super- 
conducting NbTi filaments, chosen for their extremely 
low thermal conductance. Thermal links from the di- 
lution refrigerator to the sample are made from Man- 
ganin or PtW wire and the time constant of equilibration, 
T2, is chosen to be greater than 1 hour at all tempera- 
tures. The time constant of internal thermal relaxation, 
Ti, is always much shorter than the external time con- 
stant, t i , thereby minimizing thermal gradients within 
the sample. Twisted pairs, extensive heat sinking of leads 
and 7r-filters are employed to avoid self-heating in the 
thermometer. The sample is suspended from thin nylon 
threads (~ 10 /im diameter) inside a copper radiation 
shield that is heat sunk to the mixing chamber. The ad- 
dendum, due primarily to weak link and chip resistors, 
is determined to be < 0.1% of the total heat capacity 
in any of the conditions presented here. The thermome- 
ter resistance is measured with an LR700 ac resistance 
bridge. 



IV. RESULTS 
A. Magnetic Susceptibility 

Detailed frequency scans of the magnetic susceptibility, 
x(f) = x'{f) ~ have been taken at various tem- 

peratures in between ~ 80 mK and ~ 350 mK, on three 
different stoichiometries at Ho concentrations of 1.8%, 
4.5% and 8.0%. These scans, shown in Fig. [2] exhibit 
qualitatively what is expected of glassy relaxation, in all 
three samples. Broad (full width at half maximum rang- 
ing from 1.4 to 3 decades) and largely symmetric peaks 
are seen in the absorption spectra, x"(f)i coinciding with 
an infection point in x'(/). The peak position, defined 
by /Max; drops sharply with decreasing temperature. 

The temperature ranges studied have been roughly 
chosen as those temperatures for which the frequency 
window of our measurement (from 1 mHz to around 1 
kHz) permits the observation of the peak position in 
x" ■ Surprisingly, the resulting temperature ranges for 
all three samples are quite similar. This suggests an un- 
usual lack of scaling of the relevant temperature ranges 
with concentration, which will later be explained as a 
result of dramatically changing time scales. 

In order to compare the widths of the absorption spec- 
tra, they may be superimposed by plotting a normal- 
ized susceptibiltiy, x"(/)/XMax> a g a inst frequency scaled 
by the peak frequency, ///Max- The resulting plots are 
shown in Fig.[3j A clear, and largely symmetric, broaden- 
ing of the spectra is seen with decreasing temperature. It 
can also easily be observed that the spectra become, over- 
all, wider as the concentration of Ho 3+ ions is increased 
from 1.8% up to 8.0%. Certainly, for most temperatures, 
the spectrum of an x — 0.08 sample will be broader than 
that of an x = 0.045 sample, for example. All three sam- 
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FIG. 3: (color online) Superimposed frequency scans of the 
absorption spectra (after correcting for the demagnetization 
effect) for the three different stoichiometries, obtained by nor- 
malizing x" by the maximum, XMax an d normalizing the fre- 
quency, /, by the peak position /Max- The same legends of 
Fig. [2] apply here, with high temperatures as red and low 
temperatures as blue. A clear broadening of the spectra is 
observed with decreasing temperature. Also clear, is the in- 
creasing width of the spectra (in these temperature ranges) 
with increasing concentration, x. 



pies appear to tend toward a common shape and width 
of absorption spectrum at higher temperatures. 

Dense temperature scans of the susceptibility, that 
have been performed only on the x — 0.045 sample, are 
shown in Fig.[4j The in-phase susceptibility \' is found to 
increase monotonically as the temperature is reduced un- 
til a maximum at the frequency-dependent freezing tem- 
perature T/(/), below which it drops out sharply. x'{T) 
is shown at four frequencies of measurement, with Tf(f) 
showing a comparable behavior to the inverse of the func- 
tion fuax(T). 

A commonly employed parametrization of glassy re- 
laxation, often referred to as the Mydosh parameter,^ 
is obtained through the ratio of the fractional change of 
freezing temperature to the logarithm of the change in 
measurement frequency, thus £ = ATy / (Tf A In /) where 
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FIG. 4: (color online) Temperature scans of the 
demagnetization-corrected, in-phase susceptibility, x'(T), of 
an x — 0.045 sample of LiHo I Yi_ I F4 at four different fre- 
quencies of measurement: 0.001 Hz (purple), 0.01 Hz (blue), 
0.1 Hz (green) and 1 Hz (red). Also shown, for compari- 
son, are experimental data from Jonsson et aZP (which we 
have corrected for demagnetization assuming negligible x")> 
experimental data from Reich et ali^, numerical simulation 
data from Biltmo and HeneliusS (with an arbitrary scale fac- 
tor) and a T~ ' 75 power law proposed by Ghosh et o/P^. 



Tf represents either an average or limiting value of Tf(f). 
One could equally well apply this approach to the fre- 
quency scans of Fig. [2j by changing the definition to 
£ = AT/(TAln/ Max ). The resulting values, when ap- 
plied to the very low end of our temperature (frequency) 
range, are 0.13, 0.11 and 0.11 for the 1.8%, 4.5% and 
8.0% samples, respectively. Clearly the value of this pa- 
rameter is dependent on the choice of temperature range 
employed, but serves as a rough order-of-magnitude es- 
timate. The canonical spin glasses tend to exhibit val- 
ues of £ in the range of 0.005 to 0.06 and, as a rule of 
thumb, larger values have been taken as an indication 
of superparamagnetism or a lack of a finite-temperature 
glass transition, T 9 P^This parametrization, if taken to 
be relevant, does not strongly support the idea that our 
samples of LiHo a; Yi_ 2 ;F4 are spin glasses. However, we 
argue that the parameter £ is not a measure of the exis- 
tence of a finite T g but is, rather, a measure of how close 
to T g one is able to measure. If a measurement were 
limited to frequencies in the MHz range, instead of the 
mHz range, even the canonical spin glasses would exhibit 
a large value of £ since relaxation with such high frequen- 
cies would occur very far away from T g . In other words, 
the Mydosh parameter relies heavily on the assumption 
that all systems have similar microscopic spin flip rates 
and that all measurements are performed on similar time 
scales. We will show that this is not a safe assumption 
for the present system. 
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1. Dynamical Scaling 

Perhaps a more quantitative way of determining 
whether samples of LiHo a; Yi_ 2: F4 exhibit a spin-glass 
transition or not is to perform a dynamical scaling anal- 
ysis where 

The choice of how to define t(T), experimentally, is not 
entirely clear. Perhaps the most rigorous approach would 
be to employ the limit lim w _j. ^'(^OA^x'^)- How- 
ever, this limit is difficult to reach for most tempera- 
tures, generally requiring immensely low frequencies of 
measurement!^ Several research groups have also em- 
ployed the maximum in x'{T), or the freezing tempera- 
ture Tf, as a function of f reque ncy, in order to quantify 
the dynamics of a systemPSSSl Defining a maximum in 
temperature is a slow process, however, and also restricts 
measurements to rather high temperatures in this partic- 
ular system. We have, therefore, chosen to parametrize 
the time scales of LiHo 3; Yi_ 2 ;F4 with the maximum of 
X"(f), thus / Max , or equivalently, T Max = l/(27r/ Max ). 
These results are exhibited in Fig. [5] 

Such a dynamical scaling analysis has already been 
presented for our x — 0.045 sample in Ref. 5. Below 200 
mK, a power law of the form of Eq. ([8| was successfully 
fit to the data. With no fixed parameters, the results 
of such a fit are T g = 42 ± 2 mK, zv = 7.8 ± 0.2 and 
To = 16 ± 7 s. The critical exponent, zv, so obtained, 
matches extremely well with the value determined from 
Monte Carlo simulations of a short-range 3d Ising spin 
glass modeP^ and also with values determined for several 
canonical spin glass systems JSEESl This provides strong 
evidence that there is indeed a spin-glass transition in 
the system. It is important to note that the intrinsic 
time constant of the system, tq ~ 16 s, is exceptionally 
long compared to most other spin glass systems that have 
been studied. For example, in Eu0.4Sr0.eS, a time con- 
stant of To — 2 x 10~ 7 s is measured, thus it is about 
8 orders of magnitude faster Such a r as is observed 
here is largely unprecedented, although very recently a 
time constant of To — 0.01 s has been measured in an- 
other Ising spin glass Dy K Yi_ x Ru2Si2p^ 

This extraordinarily long time constant implies that, 
on reasonable time scales of an experiment, it is impos- 
sible to approach the glass transition and maintain equi- 
librium. The lowest temperature data presented here, at 
77 mK, are only at a reduced temperature t = 0.8. Be- 
cause of the very large exponent zv, in order to obtain 
equilibrium data at t = 0.1 (10% above T g ), for example, 
one would need to perform a measurement with a time 
constant of around 32 years. This explains why the My- 
dosh parameter £ is larger here than in most spin glasses: 
quite simply we are far from the transition so the nor- 
malized slope of r(T) is still relatively small. Above 200 
mK the dynamical scaling law no longer seems to hold. 
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FIG. 5: (color online) Dynamical scaling plots showing that r 
plotted against (T — T g )/T g on a log-log scale can be well fit 
by a straight line. The increase in to with decreasing x can be 
seen in the horizontal shift of the lines. The fits obtained with 
no free parameters are shown for x — 0.045 and x = 0.080, 
where for x = 0.018, we have shown the fit obtained with 
fixed zv = 8. 



However, since 200 mK corresponds to a reduced tem- 
perature t ~ 4, this is perhaps not surprising and we are 
simply reaching the limits of critical behavior. 

In order to gain further insight into the reasons for 
such an anomalously long time constant, we have ex- 
panded our measurements to adjacent stoichiometries. 
The x = 0.080 sample in fact lends itself equally well 
to such a treatment. A successful dynamical scaling fit 
gives, again, a reasonable exponent of zv = 7.8. The 
glass temperature obtained is T g — 65 ± 3 mK, showing 
the expected increase in T g with concentration x. The 
intrinsic time constant, To ~ 0.1 s, is actually quite a bit 
smaller than that of the 4.5% sample, thus there appears 
to be an inverse correlation of To with x. 

The 1.8% sample does not present such a clear picture, 
however. A fit with no free parameters yields T g = 41 
mK, zv = 7.3 and To = 61 s. The critical exponent is 
somewhat small and, most importantly, the glass transi- 
tion, so obtained, barely differs from that of the 4.5% 
sample, with more than twice the average interaction 
strength. Fixing the critical exponent at zv — 8.0 also 
results in a perfectly adequate fit, with T g ~ 35 mK 
and t = 560 s. From a \ 2 analysis, the fit is found to 
be quite under-constrained, with perfectly acceptable fits 
obtained with a range of glass temperatures from 33 to 46 
mK. This may be a result of only being able to measure 
far from the glass temperature. 

Even though the error bar on T g (x = 0.018) is rather 
sizable, all the analyses above give surprisingly high 
glass temperatures. If one were to assume that T g 
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should scale linearly with x, we would expect T g ~ 17 
mK. A dependence not precisely linear, but not far 
from linear either, would be expected: some spin glass 
systems^ scale roughly as x 7 , the mean-field the- 
ory of LiHo a; Yi_ 2 ;F4 has a slightly higher than linear 
dependence^ and Monte Carlo simulations of Tarn and 
GingrasP suggest a roughly linear dependence on x. We 
are therefore led to speculate that there is either some 
reason that T g is unusually high for the 1.8% sample or 
that there is some additional complexity to the curve 
t(T) which obscures the true spin glass temperature. 

What is very clear is that To of the 1.8% sample is 
much longer than that of the 4.5% and 8.0% samples, 
again highlighting the inverse correlation of tq with x. 
In Section |VJ we discuss a possible explanation for this 
inverse behavior and the exceptionally slow dynamics in 
these materials. 



2. Shape of Spectra 

These materials clearly present challenging measure- 
ments, requiring very long time scales. But with these 
slow dynamics comes an advantage in that one can re- 
solve a larger part of the spectra with susceptibility ex- 
periments. For most spin glasses, the frequency window 
of a typical ac susceptibility measurement corresponds to 
temperatures quite close to the glass temperature. The 
resulting absorption spectrum, consists of a fea- 

ture that is immensely broad, covering many decades. 
One is able, for certain temperatures, to resolve the peak 
positio n but n ot say much more about the shape of that 
feature. ^ 47 * 48 ! Because the dynamics are inherently very 
slow in LiHo x Yi_2;F4, ac susceptibility measurements 
have a frequency window that corresponds to temper- 
atures further away from T g , where the absorption spec- 
trum is quite a bit narrower. This permits us to describe 
the shape of the spectra in more detail, including the 
widths and the power-law behavior of the low- and high- 
frequency tails. 

Most notably, all of the spectra, on all three sam- 
ples, appear to tend toward a common form in the high 
temperature limit. This form has a full width at half- 
maximum (FWHM) of roughly 1.4 decades. The low- 
frequency tail appears to be linear in frequency. In con- 
trast, the high-frequency tail seems to follow a shallower 
power law cx with f3 roughly 0.7 to 0.75. As the 
temperature is lowered, the high-frequency power law 
becomes less steep. At the lower end of our temper- 
ature range, the 4.5% sample, for example, exhibits a 
high-frequency tail with an exponent closer to f3 ~ 0.65. 
Although the 8.0% sample broadens much more at the 
temperatures studied here, the high-frequency tail does 
not get much shallower. Meanwhile, the 1.8% sample 
shows very little overall broadening as the temperature 
is reduced, yet has a high-frequency tail with an exponent 
close to /3 ~ 0.55. It is difficult to determine whether the 
low-frequency tail also becomes less steep or whether it 
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FIG. 6: (color online) The normalized absorption spectra of 
all three samples studied, when scaled on the frequency axis 
with the exponent 1 — T g /T are found to collapse onto a sin- 
gle curve, for much of the frequency range. The tails of the 
absorption spectra are not accommodated with this scaling 
relation. The same values of T g as were found in a dynamical 
scaling analysis are used to scale the 4.5% and 8.0% samples, 
whereas the 1.8% sample required a T g = 30 mK to obtain 
good overlap. Colors of the data points correspond to those 
of Fig. [5] and range from ~ 200 mK (red) down to ~ 80 mK 
(blue). All data have been corrected for the demagnetization 
effect. 



is simply that the spectra become very broad and we are 
unable to observe the limiting behavior on that side. 

These details make it difficult t o fit s uccessfully to the 
spectra. The Davidson-Cole formJ^lS3 for example, can 
reproduce the low and high-frequency power-law behav- 
ior, but cannot account for the widths of the spectra, al- 
though it comes quite close at higher temperatures. The 
ac susceptibility resulting from relaxation processes can 
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be described by the Debye model 

p(r)d\nT 



components of the hyperfine interaction and its contri- 
bution to the Hamiltonian will be 



XM = Xo 



1 



(9) 



where p is a distribution of relaxation times of the system. 
An absorption spectrum that is linear at low frequency, 
as is found in the Davidson-Cole form, results from a cut- 
off at long t (low frequency) in p(r). The distribution 
p must also have a long low r (high frequency) tail with 
power law f3 in order to give rise to that same exponent 
in the susceptibility at high frequencies.^ 

In looking at Fig. |3j it can be noticed that the width of 
the spectra (on a log scale) is, at least roughly, correlated 
with the distance from the transition at T g . Since the 
spectrum widths appear to tend toward a constant value 
at high temperatures, we make the ansatz that the excess 
width at lower temperatures is simply inversely propor- 
tional to the ratio T/T g . This can be tested by plot- 
ting x'VxMax against (f / f Ma x) {1 ~ Ta/T) which is shown 
in Fig. |6ja). For the 4.5% and 8.0% samples, a good col- 
lapse of the data over a large range of frequencies is ob- 
tained, using the same values of T g as were obtained from 
the dynamical scaling analysis. In fact, those values of T g 
seem to optimize the overlap of the data, lending further 
evidence that they are indeed the correct transition tem- 
peratures. The 1.8% sample is again less straightforward. 
A choice of T g = 30 mK, allows a good overlap. 

Alternatively, as shown in Fig. [6](b) , we can plot the 
full widths or half widths (at half the maximum) of the 
spectra in log of frequency. Here we define the maximum 
frequency as /Max and the left and right half maxima as 
occurring at /_ and / + respectively. Thus the full width 
at half the maximum (FWHM) is given by log 10 (/+//_). 
The half widths can be given by log 10 (/+//Max) and 
l°§io (/Max/ /- ) ■ Plotting all of these quantities as a func- 
tion of T g /T shows that the widths of the spectra do 
indeed follow a very simple behavior. While quite suc- 
cessful for the FWHM and much of the peak in x" , this 
scaling does not permit good agreement in the low and 
especially the high-frequency tails of the spectra. This 
is again illustrative of the unusual, shallow tails of the 
spectra. 



B. Specific Heat 

It has been understood for some time^ that the low- 
temperature specific heat of LiHo x Yi_ x F4 is dominated 
by the nuclear contribution, the result of a strong nu- 
clear dipole hyperfine coupling to the electronic moments. 
Note, the nuclear electric quadrupolar coupling is rela- 
tively insignificant!^. The resulting Schottky-like feature 
is quite large as a result of the number of degrees of free- 
dom afforded by the I = 7/2 nuclear moment. If the 
electronic moments are considered to be perfectly classi- 
cal or Ising-like, the moments can be written simply as 
the spin-1/2 variables S*. There are then no off-diagonal 



U H F = J2 A U m i S i 



(10) 



where A\\ = 2(J Z )A. The summations over the m» vari- 
ables in the partition function, Z, can simply be rewrit- 
ten as summations over Xi = rriiSf in order to write Z 
as the product of an electronic part Zq and an electronu- 
clear or hyperfine part, Zhf, defined as 



Z 



HF 



7/2 

N £ e 

x=-7/2 



A n x/2k B T 



(11) 



Then, since all thermodynamic quantities are derived 
from In Z = In Zq + ln Zhf, we can, for instance, take the 
specific heat to be a summation of a nuclear hyperfine or 
single-ion contribution C$i and a contribution from only 
the electronic magnetic moments, C m . 

However, the situation described above is an oversim- 
plification of the problem. In fact, transverse components 
of the nuclear hyperfine interaction do come into play 
if one performs a diagonalization of the full single-ion 
Hamiltonian including crystal field, hyperfine and Zee- 
man Hamiltonians, Hcf+Uhf+Hz, a 136 x 136 matrix 
(8 nuclear energy levels x 17 crystal field levels). The 
form for C$i used here has been obtained by perform- 
ing such a diagonalization. That said, with transverse 
hyperfine effects included, the partition function cannot 
be factored and a subtraction of the single-ion specific 
heat is no longer fully justified. Nonetheless, we will con- 
tinue to present data after such a subtraction in order 
to compare with other experimental works and to obtain 
an approximate measure of the system in the absence of 
nuclear moments. 

The measured, total specific heat of four different sto- 
ichiometrics of LiHoa;Yi_ x F4 is shown in Fig. [7£a). The 
specific heat is, as expected, clearly dominated by nuclear 
effects. The single-ion specific heat, Csi, is shown as the 
black curve. Performing a subtraction of C$i leaves us 
with broad peaks, as seen in Fig. [7|^b) . A broad peak 
in the specific heat is expected in spin glass materials as 
the corresponding critical exponent a is generally quite 
negative, in the range —2 to — 4) 57 * 69 l However, one does 
expect the maximum to be close to, often 20% above T g . 
This is clearly not the case here, with the peak position at 
~ 120 mK in all samples, independent of x and therefore 
independent of T g . 

It can also easily be noticed that the magnetic specific 
heat C m is overall much smaller for lower concentration 
samples. Since the entropy may be obtained from the 
relation S = (C/T)dT, there is obviously much less 
entropy accounted for in the x = 0.018 sample than in 
the x = 0.12 sample, for example. This suggests that 
there is significant residual entropy at the lowest temper- 
atures of our measurement. What happens below that 
temperature (~ 75 mK) decides what the true ground 
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FIG. 7: (color online) (a) Total specific heat measured on 
four samples of LiHosYi-^ at x = 0.018, 0.045, 0.080 and 
0.12. Also shown are the calculated single ion contributions. 
The dashed line assumes an Ising doublet with no crystal field 
excitations. The solid curve is a diagonalization of the crystal 
field and hyperfine Hamiltonians, thus includes some coupling 
with crystal field levels induced by the transverse components 
of Hhf- The latter is subtracted from the total specific heat 
to yield the specific heat of the magnetic moments, C m (T), 
which is plotted in (b). 



state, residual entropy Sq might be. Almost 100% of the 
expected i?m2 entropy is reproduced in the x = 0.12 
sample, and a steeper than linear extrapolation to T = 
(common in spin glasses^) is required in order to not 
exceed i?ln2. However, the other samples do seem as 
if they might possess varying levels of T — residual 
entropy. 



V. DISCUSSION 

While the dynamical measurements presented here 
have provided appreciable evidence of spin glass physics, 
several puzzles remain to be solved. Of utmost impor- 
tance is an explanation for the incredibly slow dynamics 
that are observed well above the glass transition tempera- 
ture T g . The specific heat is also rather confusing because 
of the stationary peak position as x is changed. Finally, 
the T g obtained here for the 1.8% sample seems anoma- 
lously high when compared to values for the other stoi- 
chiometries and even the Curie temperature, xTq{x = 1). 

We first attempt to understand the long time scales in 
this system, which are unprecedented in spin glass ma- 



terials. A possible explanation may come from a link 
to so-called superspin glasses!^ In such materials, con- 
sisting of randomly interacting magnetic particles, dy- 
namics have been observed that a re no ticeably slower 
than those of canonical spin glasses E2IZ1I Since the mag- 
netic units in a superspin glass system are often much 
larger (for example^! ~ 300/is) than the single-spin units 
of a spin glass, the microscopic spin flip time can be 
much longer. On the other hand, some objects defined 
as superspin glasses maintain rather fast spin flip times 
and values of t comparable to those of canonical spin 
glasses.^ Other systems of random magnetic particles, 
known as superparamagnets, show slow relaxation due 
to large anisotropy barriers, but not necessarily a finite- 
temperature freezing transition as in spin glasses!^ One 
could imagine that in LiHo 2; Yi_; E F4, there might be large 
ferromagnetically correlated regions that act as magnetic 
particles or clusters, resulting in superspin glass or super- 
paramagnetic behavior. However, one should expect such 
an effect to get stronger as one nears the ferromagnetic 
phase at higher x\ here we see a marked speeding up of 
the dynamics as x is increased. Furthermore, the mo- 
ments are evenly distributed throughout the sample in 
LiHo a; Yi_ a ;F4 and the samples are extremely high qual- 
ity single crystals with no indications of clustering.^ Most 
likely, these material should be treated as spin glasses 
rather than superspin glasses, with single Ho 3+ moments 
as the fundamental building blocks. 

A probable explanation for the long time scales in- 
volved can be found in single ion physics. Most im- 
portantly, the Ho 3+ ions in LiHo x Yi_ 2; F4, when consid- 
ered individually, have a truly Ising doublet ground state 
wherein the matrix elements (j" \J+\ \) and (t |J_| \) 
vanish. The result is that spin flips are forbidden transi- 
tions and must occur via the next excited state at ~ 11 K, 
or must occur via quantum tunneling. In the tempera- 
ture range of our measurements, at 200 mK and below, 
there is very little probability of exciting an ion into the 
1 7) state, thus we are left with almost entirely tunneling, 
which has the potential to be quite slow. A similar argu- 
ment might apply to the Ising material Dy x Y 1 _. I .Ru2Si2 
which shows a similarly long To of 0.01 sP^ 

In the absence of off-diagonal contributions to the 
Hamiltonian, this tunneling would not occur at all. How- 
ever, in a real material there are always mechanisms for 
spin relaxation. The tunneling here is primarily gov- 
erned by the off-diagonal or non-secular components of 
the dipolar interactio n 1 19 * 30 ! 
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where the coefficients Cij and dij are easily obtained from 
Eq. Q and x' is the projection of r,-j in the xy-p\ane. Es- 
sentially, there is another criterion for strong Ising char- 
acter of a material, that is that the dipolar interaction 
does not introduce a large degree of mixing with the next 
excited crystal-field state and therefore does not intro- 
duce a large degree of quantum tunneling. This may be 
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FIG. 8: (color online) (a) The time constants TMax as a function of temperature for the three samples of LiHo a; Yi_ :E F4 studied 
here. The solid lines are the dynamic scaling fits, described in the text, whereas the dotted lines are simply guides to the eye. 
It can be noticed that the time constant is inversely proportional to concentration x at temperatures above 150 mK. Below 
that point, the curves begin to cross as time scales in higher x samples begin to diverge with a power-law behavior at higher 
temperatures, (b) The energy splittings versus transverse field of the four lowest lying time-reversed pairs of electronuclear 
states, obtained by diagonalizing the single ion Hamiltonian Hsi = Hcf + T-Lhf + Tiz (solid lines). The dashed line is the 
energy splitting in the ground-state doublet with no nuclear hyperfine interaction, (c) An energy level diagram showing the 
16 electronuclear states. The horizontal distance between pairs is qualitatively representative of the strength of the tunneling 
matrix elements. At low temperatures, all but the lowest two electronuclear states are depopulated and relaxation is a result 
of tunneling between these states. 



reduced to the requirement that (f \Hd\i) *C 1. This 
requirement is satisfied in LiHo. E Y 1 _.j.F4 by virtue of the 
form of the states | f), | \) and The pyrochlore sys- 
tem Tb2Ti207, is an important counterexample. Despite 
an energy gap A at least as large as that of LiHo K Yi„ 2: F4, 
the coupling with those excited states through % D is 
much more significant, introducing quantum fluctuations 
and possibly making Tb2Ti207 an interesting quantum 
spin ice material.^ 

Even once assured that this system is strongly Ising, 
there is some tunneling that allows spin flips to occur 
and gives rise to dynamics at low temperature. Here, we 
engage in some basic calculations with the goal of under- 
standing whether the immensely slow dynamics observed 
in dilute LiHo x Yi^a;F4 can be realistically attributed to 
largely single-ion physics. We concentr ate here on the 
first and dominant Jf JJ term of Eq. (12). In order 
to greatly reduce the difficulty of the problem, we an- 
alyze a single Ho 3+ ion and incorporate that dominant 
off-diagonal term of TId as a transverse magnetic field. 
At any given time, nearby spins may be thought of as 
generating a random and fluctuating transverse field. We 
approximate the level of this transverse field in two dif- 
ferent ways. (1) We simply take the field generated by 
a nearest neighbor Ho 3+ ion, h^ n . (2) We take random 
distributions of ion positions and spin orientations, in 
other words snapshots of the system in time, calculate 
the transverse field at site i as 



(13) 



and average \hj~\ over many different random configura- 
tions to obtain an average transverse field h 1 - . 



While the importance of random fields in 
LiHo 2; Y 1 _ ;r F4 has been discussed extensively] 43 * 44 ! 
it has generally been static random fields induced by 
an applied external transverse field H±. Since we are 
discussing the system in zero transverse field, above the 
spin-glass transition, the random fields are fluctuating 
in time along with the spins. As such, they may not 
strongly affect the equilibrium properties of the material, 
but they certainly should have a powerful effect on the 
dynamics of the system. 

Essentially, the tunneling rate of an ion may be related 
to the energy splitting AE of the Ising doublet induced 
by the transverse field P3 Specifically, we can say that the 
time constant for a single spin flip can be approximated 
a s Tgp = h/ AE. Since the dipolar interactions in a lower 
concentration (x) sample are on average lower in energy, 
so are the random transverse fields seen by the magnetic 
ions. Thus we can immediately see why there is an in- 
verse correlation between x and r, as seen at the higher 
temperatures in Fig.[8|a). Taking, for instance, the r(T) 
curves for the 1.8% sample and the 8% sample, it can 
easily be seen that r is much larger for the 1.8% sam- 
ple at higher temperatures. The glass temperature T g , 
however, is not inversely correlated with x. Thus as the 
temperature is lowered and T g (x — 0.08) is approached, r 
of the 8% sample begins to diverge and there is a crossing 
of the curves. 

Qualitatively, the above description seems rather plau- 
sible, but is it quantitatively sufficient to explain the 
magnetic properties of LiHo^Yi-a^? First we consider 
the single-ion Hamiltonian, ignoring the nuclear hyper- 
fine interaction, thus keeping the crystal field and trans- 
verse field Zeeman interaction: %gf +Hz- We then ex- 
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actly diagonalize this 17x17 Hamiltonian and determine 
the energy splitting between the two lowest lying states 
as a function of transverse field, shown as the red dotted 
line in Fig. j8^b) . The transverse field induced by a near- 
est neighbor spin has magnitude h^ n — 0.057 T, thus 
leading to an energy splitting of 1.8 mK and a micro- 
scopic tunneling time of tsf — 4-2 ns. Taking the aver- 
age transverse magnetic field generated by a random oc- 
cupancy and spin configuration, gives h^- = 0.026, 0.052 
and 0.081 T for x = 0.018, 0.045 and 0.08 respectively. 
It seems clear that the tunneling rates determined in this 
way are far too fast to explain these materials' behavior. 

To go further, one must introduce the energetically 
significant nuclear hyperfine interaction. In several por- 
tions of the phase diagram of LiHo 2 ,Y 1 _ ;I ,F4, it has al- 
ready been observed or theoretically explained that the 
nuclear hyperfine interaction blocks the effec t of an ap- 
plied transverse field at low temperatures P22l Also shown 
in Figure j8^b) , are calculations performed by diagonaliz- 
ing the full single-ion Hamiltonian in a transverse field, 
thus Hcf + Hhf + Hz, a 136 x 136 matrix. The nuclear 
hyperfine interaction splits each of the Ising states into 
eight clectronuclcar energy levels, separated by roughly 
200 mK, as shown in Figure [8^c) . In the absence of mag- 
netic field, there are degenerate pairs of time-reversed 
electronuclear states, the lowest of which consists of 
| ^,—7/2) and | |, +7/2). These pairs are split in en- 
ergy by a transverse magnetic field, as shown for the four 
lowest- lying pairs in Fig. ^b). When compared to the 
calculation in the absence of HhFi on e can see that the 
lowest energy pair experience a drastically smaller energy 
splitting, A-B, hence dramatically slower tunneling. Cer- 
tainly transitions between higher energy electronuclear 
states are either as fast or faster than spin flips without 
the nuclear hyperfine coupling, but these states will be 
very heavily depopulated at temperatures below 200 mK 
or so. 

In a perturbative approach, one can describe the dy- 
namics in terms of virtual transitions. 30 However, there 
is certainly no energy splitting between | '[,—7/2) and 
| I, +7/2) to first-order in perturbation theory. Such a 
splitting requires many orders in perturbation theory as 
one climbs up one ladder of nuclear energy levels and back 
down on the other side of the diagram in Fig. |8^c). This 
is the reason that transitions are so immensely slow be- 
tween the lowest-lying states and the reason the nuclear 
hyperfine interaction so effectively blocks the effects of 
transverse fieldP^ 

This importance of single-ion physics was first sug- 
gested by AtsarkiniSI w ho proposed that the slow relax- 
ation seen in LiHo I Yi_ x F4 could be largely explained 
through depopulation of nuclear energy levels. Further 
evidence of such a link comes from work of Giraud et 
a ;jHHHl w ho have studied samples that are more dilute 
than those studied here by more than an order of magni- 
tude. In those samples, for temperatures above 100 mK, 
one is primarily observing single-ion properties with very 
little impact of the surrounding ions and therefore rela- 



tively minor effects of disorder. In such an experimental 
situation one is able to observe resonant tunneling of the 
electronic spins where nuclear energy levels are coinci- 
dent with the Zeeman energy gj/isH • J EMU] Where the 
system is off-resonance, the dynamics become very slow. 
Even at temperatures as high as 1 K significant absorp- 
tion occurs with frequencies as low as 1 Hz. Clearly the 
picture becomes progressively more complicated as one 
introduces larger and very random interaction strengths. 

The behavior of the real system ought to be somewhat 
challenging to model accurately. There are transitions 
between many different states to consider, longitudinal 
fields will also play a role, cotunneling of spins may be 
important and eventually the system must be thermally 
linked to a bath of phonons. However, this simple single- 
ion picture can hopefully give an estimate of what the 
tunneling rates will be at low temperatures. Taking the 
average internal transverse fields quoted above, we obtain 
estimated tunneling times of t$f = 20, 9.8 and 6.1 ms for 
x = 0.018, 0.045 and 0.08 respectively. This more than 6 
order of magnitude increase in time constant as a result 
of the hyperfine interaction must be the explanation for 
the exceedingly long time constants in this system. 

There remains, however, an appreciable discrepancy 
between these very simply calculated values and the mea- 
sured time constants, r . First we note that there is no 
particular reason that To obtained from a dynamical scal- 
ing analysis should be representative of a single tunnel- 
ing event. If that were the case, one would have a r 
much faster than the microscopic tunneling time as soon 
as the temperature was brought above T = 2T g , which 
is quite unphysical. Furthermore, it is evident that we 
have not yet reached the regime below which there is no 
temperature dependence to the tunneling rate. This is 
clear from the fact that To is exponentially dependent on 
x. In other words, the microscopic tunneling rate is fol- 
lowing an exponential Arrhenius law in temperature and 
since the critical region of a lower-x sample is at lower 
temperatures, To is exponentially dependent on x. Such 
estimates and qualitative conclusions would be greatly 
enhanced with precise theoretical work relating micro- 
scopic spin flip times with macroscopic critical behavior, 
calculations which have, to our knowledge, not yet been 
performed. Regardless of the exact quantitative result, 
it is clear that the nuclear hyperfine interaction, at the 
lowest temperatures, leads directly to a slowing of the 
dynamics by many orders of magnitude. 

While unprecedented in spin glasses, we note that such 
extreme tunneling times have been observed in the spin 
ice material Dy2Ti2 07 where the microscopic tunneling 
rate is found to be several ms" 1 even at 5 K or hi gher, 
even without a relevant nuclear hyperfine coupling! 78 ^ 79 ! 
Ho 2 Ti 2 C)7 has similarly long relaxation tiniest The pri- 
mary reason that spin ice materials are so slow is that 
the next excited crystal field energy levels are at close 
to 300 K (as opposed to 10 K in LiHo a; Yi_ 2 ;F4). Since 
those materials are not randomly diluted, the internal 
transverse fields may at times be canceled by symmetry 
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leading to even slower dynamics. 

Two main puzzles remain in our LiHo a; Yi_ 2: F4 results. 
(1) The 1.8% sample has an anomalously high T g (or ap- 
parent T g ), as compared to the other samples studied. 
Perhaps the actual T g is really enhanced at low x for 
some currently unknown reason. Alternatively, T g may 
simply appear higher. For the dynamical scaling picture 
to hold, the t$f must be somewhat fiat in temperature. 
Measurements on the 1.8% sample may be in a regime 
where there is still a temperature-dependent population 
of higher nuclear energy levels. The purely single-ion 
physics is already quite rich as a result of the nuclear 
hyperfine interaction and should consist of several Ar- 
rhenius laws before interactions take over at low T and 
power-law behavior begins to dominate. This rich nu- 
clear contribution to the dynamics may complicate the 
interpretation of the T max (T) curve and lead to an inac- 
curate determination of T g . 

(2) The other main puzzle remaining is the peak posi- 
tion of the specific heat curves and its lack of scaling with 
x, which contrasts with numerical treatments of the prob- 
lem.^ This again may be connected to single-ion physics 
and the very important coupling to the nuclear moments. 
It is certainly likely that the subtraction of Csi from the 
total C does not reveal the behavior that would be ex- 
pected without the hyperfine coupling. At such low tem- 
peratures, the system involves electronuclear states that 
are not separable.^ The peak seen at ~ 120 mK in all the 
samples is likely representative of that coupling. Perhaps 
at much lower temperatures, closer to T gi another peak 
could be found in C. However, as with susceptibility, 
specific heat measurements will eventually also fall out 
of equilibrium thus performing measurements at lower T 
than presented here is quite challenging and potentially 
misleading. 



VI. COMPARISON WITH OTHER GROUPS 
A. Antiglass 

In work of Reich et q/LLLL ZI an( j later, by the same re- 
search group, Ghosh et aL, | n | 121 a highly unusual state 
was observed at a concentration of x = 0.045, which 
has been referred to as the "antiglass" state, as a re- 
sult of certain properties that are quite the opposite of 
those expected for a spin glass. The phenomenology of 
the antiglass phase includes a narrowing and markedly 
increasing asymmetry of the absorption spectra as the 
temperature is reduced, a T~ a75 power law of the dc 
limit of the susceptibility and sharp peaks in the spe- 
cific heat. In contrast, we are unable to reproduce the 
above effects, and find more conventional behavior that 
is consistent with a spin glass, for x — 0.045 and two 
concentrations surrounding it. 

While the antiglass state was characterized by absorp- 
tion spectra that narrow with decreasing T and develop 
a strong asymmetry, we see clear and largely symmetric 
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FIG. 9: (color online) Comparison of absorption spectra mea- 
sured at 120 mK, showing good agreement between that mea- 
sured in this work and by Reich et al^ but anomalous results 
from Ghosh et a/P^ Note that no arbitrary scaling of the 
vertical axis has been performed here. The inset shows the 
maximum frequency of \ ' '(uj) as a function of temperature 
for the same three measurements. 



broadening down to 77 mK. Furthermore, the tempera- 
ture dependence of the peak frequencies /Max(^) between 
our results and those of Ghosh et al. are very different 
(see Fig. [9]) . The dc susceptibility reported to follow a 
y-o.75 p 0wer i aWj j s m f ac t found here to be much steeper 
with a temperature dependence closer to T~ 12 for much 
of the temperature range of our data (see Fig.|4]). Finally, 
the specific heat presented in Ref. [12] (part of which was 
presented previously in Ref. fTJ , exhibits two sharp peaks 
at 300 and 100 mK, and does not resemble our measured 
specific heat which has a single broad feature at around 
120 mK (see Ref. [3]). 

We note that results by Reich et al. and Ghosh et al., 
from the same research group and on the same samplep^ 
do not match in several ways. A large part of the data of 
Reich et al. agree quite well with our data. We observe 
similar x'(f) ~ *x"(/) spectra, with a comparably slow 
relaxation rate even well above T g . The peak frequencies, 
/Max(7 1 ) are almost identical, and the dc limit of our data 
XdciT) is quite close. Reich et al. show a slight broaden- 
ing of the absorption spectra above 150 mK, which we do 
not observe, and our specific heat results disagree. Oth- 
erwise, we obtain nearly the same results as Reich et aZP 
at a Ho concentration of x — 0.045. The behavior that 
clearly points to spin glass physics in our results is found, 
for the most part, at lower temperatures and frequencies 
than were obtained by Reich et al^ 

It is has been suggested that a possible explanation for 
the differences between the work of Ghosh et aZP^and our 
results as previously published^ could come from differ- 
ences in sample geometry. The measurements presented 
here were taken on two different sample geometries, one 
needle shaped and one less elongated. The resulting spec- 
tra are in fact different, but this is likely a result of the 
demagnetization effect. In the case of the contentious 
x = 0.045 sample, we have also measured a third geome- 



14 



try (very close to the geometry measured by Ghosh et al.) 
to verify that our method of demagnetization correction 
is successful and that there are no strong changes in the 
true material x(/) associated with the sample geometry. 
As expected, a different shape of spectrum is observed 
for this third sample geometry, but the result is precisely 
accounted for by the demagnetization correction. Fur- 
thermore, the differences between the results of Ghosh 
et al. and those presented here are far too significant to 
be explained with the effects of sample geometry, as is 
demonstrated in Fig. [10] 

When comparing our results and those of Reich et al. 
to those of Ghosh et al., we see that those of Ghosh et 
al. p3 are found to have higher characteristic frequencies 
and a shallower dc \' as a function of temperature. Since 
the characteristic frequency is a monotonically increasing 
function of T and the dc susceptibility is monotonically 
decreasing in T, it appears that the measurements of 
Ghosh et al. are essentially warmer than expected. This 
may be explained through heating from larger oscillating 
fields, a problem of heat sinking, a lack of thermal equi- 
librium or perhaps a problem capturing quantitatively 
accurate susceptibility values at lower frequencies. The 
sample studied by Ghosh and coworkers has a surface to 
volume ratio at least 10 times smaller than the samples 
studied here, possibly making it much more difficult to 
adequately heat sink to the dilution refrigerator. The 
heating might also involve a complex mechanism, for ex- 
ample a phonon bottleneck, where heat dissipated into 
the magnetic moments is not easily released via phonons, 
effectively leading to poor heat sinking at low T and 
strongly nonlinear heating effects. The modified shape 
of the absorption spectra may be a result of variable 
amounts of heat dissipation in the sample at different fre- 
quencies of measurement, owing to the strong frequency 
dependence of x" ■ This distortion of the spectra, which 
occurs predominantly on the low-frequency side, seems to 
lead to an inaccurate measure of TMax(T') and an appar- 
ent narrowing of the spectra with lower T. These ideas 
remain largely conjecture, however, as we have not been 
able to simulate or experimentally reproduce the results 
of Ghosh et al. by introducing heating to the sample. 



B. Nonlinear susceptibility 

A third group of researchers, Jonsson et al., have also 
measured the bulk susceptibility of dilute LiHo^Yi^j^ 
single crystals^ These experiments involved linearly 
sweeping the external magnetic field and observing the 
resulting magnetization with a SQUID. Then, over a cer- 
tain range of magnetic field, they perform a fit of the 
form M(H) = X\H — X3H 3 , in order to extract the lin- 
ear and nonlinear susceptibility. Their linear susceptibil- 
ity result is shown in Figure [IJ and is found to match 
well with ours down to a temperature of around 150 
mKpSO Although Jonsson et al. are sweeping the mag- 
netic field quite slowly, there is guaranteed to be some 
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FIG. 10: (color online) Raw data, not corrected for demag- 
netization. Our data, obtained at 95 mK, on two differ- 
ent sample geometries: a needle-like sample of dimensions 
0.57 x 0.77 x 7.7 mm 3 (blue diamonds) and a much less elon- 
gated sample of dimensions 0.57 x 0.77 x 1.2 mm 3 (green cir- 
cles). Shown for comparison are data from Ghosh et a?P^, 
taken at a similar temperature of 90 mK, on a sample with 
dimensions 5 x 5 x 10 mm 3 . The difference between our two 
samples is well accounted for by the demagnetization effect, 
whereas the data of Ghosh et al. are quite different, shifted to 
higher frequencies and quite assymetric. Solid lines are guides 
to the eye and not theoretical fits. 



loss of equilibrium below a certain point, which appears 
to be occuring below 150 mK or so. Because of the dif- 
ferent measurement techniques employed, it is difficult to 
directly compare our results, particularly since the mea- 
surement of Jonsson et al. does not operate at a unique 
or well defined frequency. We can get a very approximate 
idea of what frequency of measurement is being employed 
through / ~ (dH/dt)/ AH, where AH is the range of 
fields over which the data are analyzed and (dH/dt) is 
the sweep rate employed. Given the numbers quoted in 
Ref.|H that frequency is, at the lowest, ~ 40 mHz. While 
indeed quite a low frequency by most standards, it can 
be seen from Figure |4j that one should see a freezing 
temperature of around 110 mK. 

Thus, our results do not seem to be in disagreement 
with the data of Jonsson et al., but rather disprove the 
idea that their measurements are in the dc limit. From 
the nonlinear susceptibility, X3 > they have found a rather 
broad peak, and have concluded that that feature is not 
sharp enough to represent a spin-glass transition. They 
have argued, on the same basis, that there is not a spin- 
glass transition in an x — 0.167 sample either, which has 
opened up a debate regarding a part of the phase diagram 
previously thought to be well understood!^ It is true 
that the appearance of a peak in X3j n ot accompanied by 
a scaling analysis, is insufficient evidence of a spin-glass 
transition. 

However, based on our dynamical measurements, we 
believe that the dynamics are rounding this peak and 
obscuring the glass transition. While the critical expo- 
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nent for X3 is 7 — 3 in most spin glass systems, indeed 
leading to a sharp divergence, the dynamical critical ex- 
ponent is much larger: zv ~ 8. Thus the dynamical 
component will always "win" at some point and eventu- 
ally round the peak in X3- How far away from the true 
T g this occurs is dependent on To and the frequency of 
measurement employed. Thus for materials with a small 
To, something close to a divergence in \3 may easily be 
obtained, but for systems with larger t , closer attention 
must be paid to the dynamics. A measurement of X3 
at low x in the LiHo a: Y 1 _ a ;F4 series, could be a valuable 
future measurement and a definitive proof of spin glass 
physics, but caution is warranted, as the extreme values 
of To should make it extremely difficult to obtain a sharp 
peak in % 3 and a careful dynamical scaling analysis would 
likely be necessary. 



/iSR experiments - Recently, Rodriguez et 
have performed /iSR experiments on the same sam- 
ples of LiHo 2; Y 1 _ :l; F4 that were used for the work pre- 
sented here. Although they do not find any significant 
change in behavior between higher concentration sam- 
ples (x = 0.12 and x = 0.25) and concentrations around 
x = 0.045 that would be suggestive of an antiglass state, 
neither do they obtain results that are typical of spin 
glass physics. In fact, the /iSR experiments show a 
temperature-independent fluctuation rate below at least 
200 mK that is faster than the lower limit of the fre- 
quency window of fiSK measurements. This is reminis- 
cent of the persistent spin dynamics (PSDs) that are ob- 
served in many geometrically frustrated magnets, includ- 
ing the pyrochlore spin glass Tb2Mo2 07)S21 



It is quite difficult to reconcile the idea of persis- 
tent spin dynamics exhibiting fluctuation rates ranging 
from 0.75 to 20 MHz with the sharply temperature- 
dependent dynamics reaching frequencies of less than 
1 x 10~ 3 Hz, seen with susceptibility measurements in 
the same temperature range. In several other magnetic 
systems, the microscopic fluctuation rate as determined 
with /iSR has been found to be much faster than the 
characteristic freque ncy o btained from bulk ac suscepti- 
bility measurements E3HS31 Included is the dipolar spin ice 
material Dy 2 Ti 2 07, in which the Dy 3+ ions also carry 
Ising moments. There, a qualitatively similar temper- 
ature dependence of the dynamics is observed in both 
^SR and in ac susceptibility, but there is an overall mis- 
match of time scales, many orders of magnitude in size.^1 
In LiHo a; Yi_ 2 ;F4, we do see some commonalities between 
the /LtSR and x"(f) data: there is a significant drop in 
the low-temperature fluctuation rate, Uq, as x is lowered, 
coincident with our observation of an increase in t . Ro- 
driguez et alWl also draw attention to the importance of 
the nuclear hyperfine interaction in this system with the 
observation that the plateau of the fluctuation rate has 
an onset near 300 mK, not so distant from the 200 mK 
energy scale associated with Hhf- 



C. Theory 

With some exceptions, our results are very close to 
theoretical work on the ideal dilute dipolar Ising model. 
The glass temperatures obtained here for x = 0.045 and 
x = 0.08 are well matched to those determined in Monte 
Carlo work of Tarn and Gingrasp^ as can be seen in 
Fig. 11 'b). Although the theoretical work determines T g 
for the concentrations x = 0.0625 and x — 0.125, we can 
easily interpolate between those points and find adequate 
agreement between experiment and theory. The 1.8% 
sample shows an anomalously high T g when compared 
to theory, however. In section V, we have suggested two 
possible nuclear hyperfine effects that may explain this 
discrepancy. Biltmo and HeneliusP have also performed 
Monte Carlo simulations on the ideal model, notably with 
x = 0.045, and have computed several physical proper- 
ties, including the dc susceptibility and specific heat. The 
dc limit of our ac susceptibility measurements are in fact 
very much in line with those simulations, as shown in 
Fig. [4] The specific heat, however, is found in theory 
to exhibit a broad peak that shifts in temperature with 
x. The measured specific heat exhibits a similarly broad 
peak, but one that does not change with x. Thus this 
discrepancy between experiment and theory remains an 
important question to issue to resolve. 



VII. CONCLUSIONS 

To summarize, we have shown that a dynamical scal- 
ing analysis of ac susceptibility measurements taken on 
low- a; samples of LiHo x Yi_ x F4 gives strong evidence of 
finite-temperature spin-glass transitions. Furthermore, 
these glass transition temperatures, T g , are well matched 
to theory^ as can be seen in the phase diagrams plotted 
in Fig. [IT] There remains some ambiguity regarding the 
glass temperature T g of the 1.8% sample and the specific 
heat of all the samples that remains to be fully under- 
stood. We have tentatively proposed that these anoma- 
lies in our understanding of this system may be the result 
of the complicated role of the nuclear hyperfine interac- 
tion on thermodynamic and dynamic quantities. Obser- 
vation of the antiglass phasep ^ 1 * 12 * 17 ! conflicts with the 
behavior that we have deduced with our ac susceptibility 
and specific heat experiments. While tentative explana- 
tions for these discrepancies have been suggested, a final 
and quantitative explanation remains to be discovered. 
Given recent theore tical conclusions ™ inconsistencies in 
the antiglass picturej^^and the clear appearance of very 
slow dynamics in ac susceptibility with conventional scal- 
ing exponents, zv, we feel, nonetheless, that there should 
remain little doubt that LiHo x Yi_ x F4 does indeed freeze 
as a spin glass at low x. 

Most interestingly, our results show compelling evi- 
dence for the very important effect of nuclear hyperfine 
interactions on the dynamics of the system. Typically 
nuclear moments are coupled weakly and have very little 
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FIG. 11: (color online) Current understanding of the phase diagram of LiHo a; Yi_ I F4 as a function of concentration x and 
temperature T. Ferromagnetic transition temperatures are obtained at x = 1 by Bitko et al. 2 , at x — 0.44 and x = 0.67 by 
Reich et al. 1 , at x — 0.45 by Quilliam et al. and at x = 0.30 by Kjaer et alW^. Spin glass transition temperatures are determined 
experimentally in this work (blue symbols) and by Anconna-Torres et alW (purple triangles) and with Monte Carlo by Tam and 
Gingras^ (green squares). A rather large error bar is shown with the poorly defined data point at x — 0.018. The solid lines 
are guides to the eye rather than theoretical models. In part (a), a region of reentrance has been drawn, as is often observed 
in other such systems.^ However, this occurrence has not yet been observed in this system. 



impact on electronic magnetism, although they can rou- 
tinely serve as a means to study the electronic physics of a 
system through NMR or Mossbauer experiments, for ex- 
ample. This system represents a rare case where the nu- 
clear moments are hugely influential, in blocking the ef- 
fects of internal, fluctuating transverse fields hence appre- 
ciably slowing the dynamics of the system. The impor- 
tance of the nuclear hyperfine interaction has also been 
observed or theoretically proposed under the influence of 
an external transverse magnetic field, when stu dying the 
effects of quantum fluctuations on the system P 16 ! 20 ! 43 ! 
These unusually slow dynamics have likely caused some 
of the confusion re gardi ng the existence of a spin glass 
state at x = 0.167,^^ where experiments with an ar- 
bitrary time scale were used to argue for or against a 
finite T g . The results obtained here show the importance 
of exploring scaling relations and a careful extrapolation 
to the dc limit, particularly when studying glassy sys- 
tems. Precise calculations of the tunneling rates of the 
Ho 3+ moments are beyond the scope of this work, but 



it is hoped that such an exercise may be pursued in the 
future. Such calculations may also have important im- 
plications for the dynamics of other rare-earth Ising sys- 
tems, including those of the recently discovered m agneti c 
monopole excitations in dipolar spin ice materials ! 55 * 78 ^ 
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